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Abstract: We analyze nonleptonic B s decays to a charmonium state and a light meson, 
induced by the b — > ccs transition, which are useful to access the B s -B s mixing phase 
jS a . We use generalized factorization and SU(3)f symmetry to relate such modes to corre- 
spondent B decay channels. We discuss the feasibility of the measurements in the various 
channels, stressing the importance of comparing different determinations of /3 S in view of 
the hints of new physics effects (NP) recently emerged in the B s sector. Finally, adopting 
a general parametrization of NP contributions to the decay amplitudes, we discuss how to 
experimentally constrain new physics parameters. 
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1. Introduction 



The detailed analysis of CP violation in particle physics is a powerful tool to test the 
standard model (SM) of elementary interactions and enveil the effects of new interactions. 
The fundamental role in the SM description of CP violation is played by the Cabibbo- 
Kobayashi-Maskawa (CKM) mixing matrix, which is unitary and implies CP violation if 
it is complex. The constraints stemming from unitarity can be represented as triangles, 
the lengths of whose sides are the moduli of products of CKM elements, while the angles 
represent relative phases between them. The most studied bd unitarity triangle is defined 
by the relation V U( [V* b + V c dV* b + V t dV t l = 0, and has been probed mainly through the ex- 
tensive analysis of Bd phenomenology. As a result, the CKM parameters in the Wolfenstein 
parameterization have been fixed with small errors through the measurement of the sides 
and the angles of this triangle [El]. The next, already ongoing, effort is to look at processes 
in which to test the SM requires a greater experimental and theoretical precision. The 
B s sector is suitable for such a purpose. As in the Bd case, one of the CKM unitarity 
constraints involves matrix elements related to B s decays: V us V* b + V cs V* b + Vt s V t * b = 0, 
and one of the angles of the corresponding bs triangle is the phase of the B s -B s mixing: 
(5 S = Arg \-¥r& ■ In SM (3 S is expected to be tiny: (3 S ~ 0.017 rad. 

L cs cb J 

B s is produced at the B factories running at the peak of T(5S) and in hadron collisions. 
In particular, the experiments CDF and DO at the Tevatron have obtained a number of 
remarkable results, such as the measurement of the mixing parameters: The mass difference 
of the two B s mass eigenstates has been fixed to Am s = 17.77 ± O.lO(stat) ± 0.07(syst) 
ps _1 [Q] , while the value of their width difference Ar s depends on the constraints on 
f3 s adopted in the experimental analysis Q; noticeably, Ar s is not small as in the B& 
case. Furthermore, these Collaborations have provided us with results which seem to 
signal new physics (NP) effects. The first one concerns the phase (5 S , extracted from the 
angular analysis of the time-dependent differential decay width in the process B s — > J/tpcfi. 
The study is rather involved: an angular analysis is needed to disentangle the CP-even 
and CP-odd components, required since the final state of two vector mesons is not a 
CP eigenstate. Moreover, the measurement can be carried out either considering flavour 
tagged or untagged decays. Another issue concerns the use or not of assumptions on 
the strong phases among the different helicity amplitudes in the considered process: this 
assumption has been once adopted by DO Collaboration in one study Q. Different results 
have been obtained from the different analyses [||, |5], and, averaging them, the Heavy 
Flavour Averaging Group has provided a value of f3 s consistent with SM only at 2.2 a 
level: = -2/3, = — 0.77±g;|f or ^l 1 ^ = -2{3 S = — 2.36drg : || f§. A new measurement 

announced by CDF: /3 S E [0.0, 0.5] U [1.1, 1.5] (at 68% CL) §, if confirmed, would reconcile 
the SM prediction with experiment. 

Another signal of a possible inadequacy of the SM is the measurement of an anomalous 
like-sign dimuon charge asymmetry of semileptonic 6-hadron decay, reported by the DO 
Collaboration || (updating a previous measurement ||): 

, 7V++ - Nr~ 

At = — ^-p h — = -(9.57 ± 2.51 ± 1.46) x 10 -3 . (1.1) 

si N ++ N v I K J 
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Hence, there is a large excess of negatively charged muons over positively charged ones 
which would have been generated by the oscillation of one neutral b meson into the other, 
at odds with the SM expectation A b sl = (-0.310lg;g||) x 1(T 3 a result which might 
imply a NP effect in the oscillation. 

In this complex scenario it is worth to further analyze the B s sector, trying to identify 
and reduce the uncertainties affecting the theoretical predictions, with the aim of improving 
the measurement of /3 S , overcoming the difficulties in B s — > J/tfj<f). Notice that this channel 
is considered a golden mode, since it is induced by the b — > ccs transition in which, in SM, 
the only weak phase involved is that of the mixing, so that the indirect CP asymmetry 
would be proportional to sin(2/3 s ), much in the same way as — > J /ipKs has provided 
a determination of the angle /?. The feasibility in the reconstruction of the products of 
the subsequent decays J/tp — > fj, + fi~, <j) — > K + K~ makes this channel also experimentally 
appealing. 

There are other modes that can be used to access (3 S , namely B s — > M( cg ) + L, where 
M( cg ) is a charmonium state J/tp, ^>(2S), r] c , r/ c (25), XcO, Xci, Xc2, h c and L is a light scalar, 
pseudoscalar or vector meson, /o(980), ry, r( and 4>. Each of these channels presents specific 
features and advantages/difficulties which we want to discuss here. Standing the general 
theoretical difficulty in the calculation of nonleptonic decay amplitudes, in the next Section 
we discuss approaches to afford the problem, and exploit the generalized factorization to 
calculate the branching fractions in the SM. In this way the most suitable processes to 
determine (3 S can be identified. In Section ||] we also consider the possible impact of new 
physics in these modes, and discuss how to exploit experimental data to constrain NP 
parameters. Conclusions are presented in the last Section. 

2. B s — > M c5 L decays 

In SM, the effective hamiltonian governing the nonleptonic decays induced by the b — > ccs 
transition reads as [II]: 



(-i (- 10,77,83 . 

n cS = -||y c6 K*[Ci(M)o 1 (/x) + a 2 (/i)o 2 ( M )] -v tb v*[ ^ 0^)0^)] j, (2.1) 

where Gp is the Fermi constant, the Cj are Wilson coefficients, and Oi are 

• current-current (tree) operators 

Ox = cr M (l - lb )b vf(l - 75 )c, 2 = C7 M (1 - 7 s)c aY*(l " 7s)b (2-2) 

• QCD penguin operators 

O3 = (s a b a ) V -A ^2(%<lp)v-A, O4 = (s/3&a)y-A ^{q' a q'p)v -A, 

q> q' 

5 = (s a b a ) V - A ^(cfpq'^v+A, 6 = (s /3 b a ) V -.A^2@a<lp)v+A, (2-3) 
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electro-weak penguin operators 

3 3 

q> q> 

3 3 

°9 = -x(s a b a )v-A^2e q '(q / pq^)v-A, Oio = -(sfib a )v-A ^2 ^'^'aQ^v-A, (2.4) 



• magnetic moment operators 

7j = -^s a a^ u (m s P L + m b P R )b a F fl „, 

Os g = -^fs a a^{m s P L + m b P R )T^b p G%. (2.5) 

a and j3 are color indices, and q' are the active q' = (u,d, s,c, b) quark fields at the scale m b 
with charge e g >. The right (left) handed current is defined as {q' ol q'p)v±A = Qalvi^ ^Ts)^; 
with projection operators Pr ; l = lj ^ 75 • Assuming CKM unitarity and neglecting the tiny 
product V ub V* s , the relation holds: V t bV t * s = —V cb V* s . 

The hamiltonian ( |2.1| ) induces the decays of SU (3) f related states, namely the decays 
of Bj, B~ and B s ; we consider the general case of the decay of a B a meson (a = u,d,s 
being the light flavour index). The simplest approach to compute the matrix element of the 
hamiltonian ( |2.1| ) between given initial and final hadronic states is the naive factorization 
approach. In such an approach, neglecting the magnetic moment operators in (2.5), the 
B a — >• M cE L amplitude reads: 

A{B a -> M C - C L) = ^=V cb V c * s a e 2 ff (v) (M (c5) | c^(l - 75 )c |0) (L\ s lfl (l - l5 )b \B a ) , (2.6) 

C\ C4 3 / C10 \ 

where a e 2 ff (n) = o 2 (/i) + a 3 (n) + a 5 (/i) and a 2 = C 2 + j^-, a 3 = C 3 + — + -e c ( Cg + -j^-J 

and 05 = 6*5 + — + —e c ( C7 + — 1 . However, naive factorization predictions are not able 

to reproduce several branching ratios for which experimental data are available. Among 
these there are B^ decays induced by the transition b — > ccs: some of these modes, which 
are of interest for the present analysis, are listed in Table |l| together with the experimental 
branching fractions. 

Several modifications of the naive factorization ansatz have been proposed. One pos- 
sibility is to consider the Wilson coefficients as effective parameters to be determined from 



experiment [13]. In principle, this implies that such coefficients are channel-dependent. 
However, some channels could be related, namely invoking flavour symmetries, so that 
universal values for the coefficients can be assumed within a certain class of modes. In our 
case, this generalized factorization approach consists in considering the quantity i n 



(|2.q ) as a process-dependent parameter to be fixed from experiment. In particular, on the 
basis of SU (3) p symmetry, B q (B u or B^) decays can be related to analogous B s decays 
induced by the same b —> ccs transition, so that experimental data concerning B q modes 
provide predictions for B s related ones. Also this method presents some drawbacks, for 
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/i c 




B~ 


< 2.1 


3.6 ±0.9 


< U.12 






B° 


1.70 ± 0.40 


2.0 ±0.6 


< 0.36(0.66 ±0.19) 


(< 2.2) 



Table 1: Experimental results for the branching fractions B(B — > MccK^) (xlO 4 ) |l2) ; results in 
parentheses are from ||. 



example the issue of rescattering in the final state and of the strong phases in the various 
amplitudes cannot be faced [jnj]. Nevertheless, it is useful from a phenomenological point 
of view, at least to understand the size of nonleptonic branching ratios. 

A different procedure to analyze nonleptonic decays is the hard-scattering approach, 
based on the assumption of the dominance of hard gluon exchange and of the suppression 
of soft mechanisms due to low energy gluon exchanges. In this approach a nonleptonic 
amplitude is expressed as a convolution of a hard kernel, computed in perturbation theory, 
with the light-cone wave functions of the hadrons involved in the decay. In this so-called 
perturbative QCD approach (pQCD) the suppression of the soft term is achieved by suitable 
Sudakov suppression factors, but the uncertainty in the wave functions limits the accuracy 
of the predictions [15]. 

A systematic improvement of naive factorization is QCD factorization (BBNS) |l6| . 
In this approach, a factorization formula is written for a nonleptonic B a — > M%M% decay 
amplitude (Mi denotes the meson picking up the B a spectator quark), valid in the heavy 
quark limit (i.e. up to h.QCD/fnb corrections). This formula reproduces the naive factor- 
ization result at leading order in a s and Aqcd/ 771 ^ however, it cannot be applied when 
the meson that does not pick up the B a spectator quark is heavy. The knowledge of the 
meson wave functions is required and represents therefore a limiting factor. 

A particular case is represented by the decays B — > M £c L considered here. Since 
the charmonium state is a heavy meson, the BBNS factorization formula does not hold. 
However, it has been pointed out that, being a charmonium meson a state with small 
transverse extension, one can still adopt the factorization formula. However a problem 
arises going beyond the leading twist for the wave functions, since the factorization formula 
contains convolution integrals of such wave functions, and higher twist wave functions do 
not vanish in the end point, developing divergences. Getting rid of such divergences requires 
the introduction of a cutoff, a parameter to be fixed from experiment. 

Two body nonleptonic B decays have also been analyzed in a modified formulation 
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of light-cone QCD sum rules originally proposed in [17] to calculate the B — > tttt matrix 
element, finding results in agreement with QCD factorization. Applying this approach to 
B to charmonium decays, one finds that nonfactorizable contributions are important, but 
that their inclusion does not allow to reproduce experimental data for B — > J/tpK |lq] . 

Hence, no satisfactory treatment of nonleptonic B to charmonium decays exists at 
present, each method having its own advantages/drawbacks. Standing the phenomenologi- 
cal importance of these modes, we afford a study based on generalized factorization, aiming 
at establishing at least the sizes of the branching ratios of these modes and their role for a 
measurement of (3 S . 

To apply eq. ( |2.6D to the modes we are analyzing, we need the following hadronic 
quantities: 

• charmonium decay constants: 

(Vc{q) |c7 M 75C |0) = -if Vc qi, , 
(J/tpfae^&y^c |0) = /^m^e* , 
(Xci(g,e)|c7 jU 7 5 c |0) = / Xcl m Xcl e* , (2.7) 

(e(A) polarization vector); for Xc0,2 and h c one has (Xco(q)\ cry^c |0) = {Xc2(q, e)| e7 M c |0) 
(h c {q,e)\c^{l-j 5 )c |0) =0; 

• B a — >■ L form factors, with L a pseudoscalar (P) or a scalar (S) meson: 



(P(SW)\s^(K)b\B a (p)) = F^q 2 ) 



2 2 

m n — mi 



(p+pT- Ba 2 P(5 V 



2 2 

+ F (g 2 ) ~ 2 mp(S V , (2.8) 



B a L form factors, with L a vector (1/) meson: 

(L(p',e)\s^(l- l5 )b\B a (p)) = 2V{ f ] W ^ 

-i [e^(m Ba +m i )^ 1 (g 2 )-(e*.g)(p + p')M 



mB + mi 

g )=^(A 3 (^)-^o(? 2 ))^] , ( 2 -9) 



with A 3 (q 2 ) = mBa +mL A 1 (q 2 ) - ^ — A 2 (g 2 ). The first two equations in M 

2m l 2m l — 

also hold for rj c (2S) and ip(2S), respectively, while the vanishing of the matrix elements 

(Xc0,2(?)| C7 M c |0) and (h c (q,e)\ cry^(l - 7 5 )c |0) implies that the B a -)■ Xc0,2 ^ and B a -> 

/i c L amplitudes vanish in the factorization approximation. 

By the factorization ansatz one has expressions for the various decay widths. Moreover, 

for decays in two J = 1 mesons, also the polarization fractions can be computed, namely 

fi, the fraction of the decay width when both the final mesons are longitudinally polarized 

. The results are the following: 

The two J = 1 mesons in the final state have the same helicity since the decaying B a is spinless. 
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modes where M cS is either a J PC 
J PC = 1 ++ P-wave Xci meson: 



1 charmonium state (J/ip or V'(25')), or a 
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+ 8m Mc _m v 
1 



(m Ba + my) 2 



}■ 



(2.11) 



G||^| 2 (a^) 2 / 2 A 1 / 2 ,,, 



ffl 



r(5 a -> m cc -v) 

{(m Ba + m v )[Af^ v 

\ / 2 2 2 \ 



lGvrmn 

t>a 



8my 



m 



2 2 > 

m M C ci m V i 



,)]«. - 



(m Ba + my) ) 

modes with a J PC = h P cg (rj c or rj c (2S)) charmonium state: 



(2.12) 



T(B a -> P cc -L) 



|2 („eff\2 f2 
■ csl 

32vrm| 



jpBa^rL 

r 



(m Pce ) 



X (™B. 



r(s a -> p cc -v) 
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327rm| a 



A, 
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P {m P - 



x A 3/2 {m\ a ,m 2 Pct ,m\) 



(2.13) 



(2.14) 



A is the triangular function, A(a, 6, c) = (a — 6 — c) 2 — 46c. Eqs. ( |2.10 ) and ( 2.13| ) apply to 
both the cases in which the light meson is pseudoscalar or scalar, L = P, S. 

Using in these expressions the coefficients ai{fx) computed in renormalization group 
improved perturbation theory, the experimental data are badly reproduced: the b —> ccs 
induced modes under scrutiny are colour suppressed, and the predictions of naive factoriza- 
tion undershoot the data. The most striking discrepancy is for the modes with XcO in the 
final state, which have a sizeable rate despite their amplitude vanishes in the factorization 
approach. Our strategy is to exploit the data in Table [l] to determine an effective parameter 
a 2 (generally channel-dependent) and, assuming SU(3)p symmetry, to use these values 
to predict the flavour related B s decays. Since the results depend on the form factors, to 
estimate this hadronic uncertainty we use two sets of form factors factors computed by 
variants of the QCD sum rule method [19], the set in JHJ obtained using sum rules based 
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on the short-distance expansion, and the set in |^] based on the light-cone expansion. In 
the case of B s — > (j) and B s — > /o(980) we use form factors determined by light-cone sum 
rules H|, Hi. 

The numerical inputs V cb = 0.0412 ± 0.0011, V cs = 1.04 ± 0.06, t{B") = (1.638 ± 
0.011) ps, t(B°) = (1.530 ± 0.009) ps, t(B s ) = (1.470±2;^f) ps, together with the val- 



ues of the meson masses, are taken from the Particle Data Group [12]. Moreover, from 



J/ip(if)(2S)) -> e+e- g we obtain: / J/v , = (416.3 ± 5.3) MeV and /^ (2S) = (296.1 ± 2.5) 
MeV. The decay constant of ij c comes from rj c — > 2j: f Vc = (380.0 ±87.1) MeV, while only 
the upper bound f Vc (2S) < 438.8 MeV is known. In the heavy quark limit, the pseudoscalar 
and vector charmonia are collected in a doublet of states with degenerate masses and same 
decay constants. Therefore, f Vc (2S) can De obtained using the relation 

f Vo (2S) = t^U(2S) = (270.3 ± 62.0) MeV; (2.15) 

JJ/ip 

symmetry breaking terms, coming from removing the meson degeneracy, are expected to 
cancel in the ratio. Since the constant f Xcl is not known, for the modes involving Xd we 
determine the product a| fx c i f rom data. SU(3)f symmetry allows to relate B s decays 
to those listed in Table [|: data on B — > M c5 K allow us to predict B s — > M C cT]('\ while 
information on B — > M c5 K* is used to predict B s — > M cS <j). As for B s —> M cE /o(980), they 
are obtained using the effective 1 0,2 1 determined from B — > M c5 K. The B s — > rf* form 
factors are related to the analogous B — >• K form factors: for a generic form factor F we 
have F Bs ^ v = — sm0F B ^ K and F Bs ~* ri ' = cos 9F B ^ K where 6 is the mixing angle in the 



flavor basis |24 



r] = r] q cos 9 — rj s sin 8, 

r\ = r] q s'm8 + T] s cos 9, (2-16) 

with rj q = (uu + dd)/\/2 and r] s = ss. The mixing angle between T] q and r] s can be fixed 
to the value measured by the KLOE Collaboration: 6 = (41.5 ± 0.3 s t a t ± 0-7 sys t ± 0.6^)° 
1 25], which agrees with the outcome of a QCD sum rule analysis of the radiative (f) — > 
modes [26]. In we include the uncertainty on the form factors at q 2 = and on 



mode 


UCDSSl 
\ a 2 


\n BZ \ 
\ a 2 1 


mode 


|„CDSS| 
\ a 2 


\n BZ \ 
\ a 2 1 




0.40 ±0.007 


0.26 ±0.005 




0.36 ±0.03 


0.25 ±0.02 


il>(2S)r)tf) 


0.50 ±0.02 


0.31 ±0.01 


Vc(2S) V ( V >) 


0.31 ±0.08 


0.21 ±0.06 


mode 


UCDSSf 1 
\ a 2 JXcil 


\a BZ f 1 
\ u 2 JXcil 


mode 


|„CDSS f 1 
\ u 2 JXcil 


\n BZ f 1 

\ a 2 JXA\ 




0.122 ±0.006 


0.076 ± 0.004 


Xcl fo 


0.122 ±0.016 


0.076 ± 0.010 


Xcl<t> 




0.0345 ± 0.006 









Table 2: Effective Wilson coefficients and combination a\ / Xol (in GeV) appearing in the 
decay amplitudes of the various B s — > M ci L modes, obtained from the SU(3)f related B decay 
modes and using the form factors in Ref.[|§ (CDSS) and Ref.|2^] (BZ). 



the experimental branching ratios in Table |T[ In the case of the transitions involving rj or 
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rf , the uncertainty on the form factors is not included since, on the basis of SU(3)f, the 
dependence on the form factors cancels when B s — > M c5 r]^ branching ratio is related to 
B —> M cc -K, leaving only a dependence on the 77-7/ mixing angle. The resulting values of 
0,2 are collected in Table ^, and the predictions for B s branching ratios in Tables ^ and 
In Table ^ also the available experimental data are included, with a general agreement 
with the predictions. Several remarks are in order. The values of a 2 derived from the 



mode B (CDSS) B (BZ) 



Exp. 



mode B (CDSS) B (BZ) 



J/ipr] 


4.3 ±0.2 


4.2 ±0.2 


3.32 ± 1.02 


llcV 


4.0 ±0.7 


3.9 ±0.6 


J/^j rf 


4.4 ±0.2 


4.3 ±0.2 


3.1 ± 1.39 


Vc rf 


4.6 ±0.8 


4.5 ±0.7 


^(25) r? 


2.9 ±0.2 


3.0 ±0.2 




Vc(2S) V 


1.5 ±0.8 


1.4 ±0.7 


1>(2S)rf 


2.4 ±0.2 


2.5 ±0.2 




Vc(2S) V ' 


1.6 ±0.9 


1.5 ±0.8 


J/i>4> 




16.7 ±5.7 


13 ±4 






15.0 ±7.8 


^(25)0 




8.3 ±2.7 


6.8 ±3.0 








XdV 


2.0 ±0.2 


2.0 ±0.2 




Xcl fo 


1.88 ±0.77 


0.73 ±0.30 


Xd rf 


1.9 ±0.2 


1.8 ±0.2 




Xd 4> 




3.3 ± 1.3 


J/ipfa 


4.7 ± 1.9 


2.0 ±0.8 


< 3.26 


Vcfo 


4.1 ± 1.7 


2.0 ±0.9 


^(25) f 


2.3 ±0.9 


0.89 ±0.36 




Vc(2S) f 


0.58 ±0.38 


1.3 ±0.8 



Table 3: Branching ratios (xlO 4 ) of the decays B s 
and in 



M cE L using the form factors in 



20| (CDSS) 

(BZ). The experimental results are taken from PDG 12 , except for B s — > J/ipr)(r]') 
measured by Belle Collaboration, with the errors combined in quadrature [^7|. For B s — > M C c/o, 
the effective coefficient 02 obtained from the B — > K mode and the form factors CDSS and BZ are 
used. The bound for B s — > J/ijj fo is due to the Belle Collaboration (2?| . 



form factors in Ref. |20| are larger than the ones derived from Ref. [21]; they also turn 
out to be channel-dependent. Their range (0.2 — 0.3) or (0.3 — 0.5) is larger than the one 
obtained in the QCDF and pQCD approaches which undershoot the data. As appears from 
Tables |3] and |], all the modes have sizeable branching fractions, so that they are promising 
candidates for measurements of /3 S . The modes involving 77, rf , /o present, with respect to 
B s —7- J/ip4>, the advantage that the final state is a CP eigenstate, not requiring any angular 
analysis. However, the channels with r/ and rf could be useful only after a number of events 
will be accumulated, since at least two photons are required for the reconstruction. 



mode 


B 


mode 


B 


mode 


B 


XcoV 


0.85 ±0.13 


Xc2V 


< 0.17 


h c rj 


< 0.23 


Xco rf 


0.87 ±0.13 


Xc2 rf 


< 0.17 


h c rf 


< 0.23 


XcO fo 


1.15 ±0.17 


Xc2 fo 


< 0.29 


hcfo 


< 0.30 


Xco4> 


1.59 ±0.38 


Xc2<l> 


< 0.10(0.62 ±0.17) 


h c 4> 


(< 1-9) 



Table 4: Branching ratios (xlO 4 ) of B s decays into p-wave charmonia. 
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As discussed in [29, 23, [3Cj, the mode B s — > J/tpfo has appealing features since, 
compared with the r\ and 7/ , the fo can be easily reconstructed in the ir + ir~ final state, 
which occurs with a large rate: B(/o — > 7r + 7r~) = (50lg)% [plf| , so that this channel could 
likely be accessed 2 . The Belle Collaboration has provided the measurement [28|: 

B(B S -> J/^/o) x B(/ -> tt+tT) < 1.63 x 10" 4 (2.17) 

which is marginally compatible with our prediction. 

The B s — > r/ c L are also predicted in Table ||. Although also these channels have sizable 
branching fractions, they present the drawback of the difficult reconstruction of the r] c . 

Let us now consider B s decays to p-wave charmonia. We have stressed that, among 
these decays, the only one with non vanishing amplitude in the factorization assumption 
is that with Xci m the final state. In the other cases, i.e. for modes involving Xc0,2 and 
h c collected in Table ^j, the results are obtained determining the decay amplitudes from 
the B decay data by making use of the SU (3) p symmetry. In this case, the differences 
between the B and B s decays arise from the phase space and lifetimes of the heavy mesons. 
As for the mechanism inducing such processes, one possibility, put forward in p2]| , is that 
rescattering can be responsible of the observed branching fractions. Among these channels, 
B s —> XcO^ is °f prime interest and very promising for both hadron colliders and B factories. 
Even though B(B S — > Xco4>) 1S one order of magnitude smaller than B(B S — > J/tp(p) it 
has appealing features, in particular as far as the potential of the LHCb experiment is 
concerned. Since the final state consists of six charged hadrons, the particle identification 
information from the RICH detectors could suppress the background. Furthermore, the 
vertex detector might be particularly efficient for these channels. 

Considering finally the polarization fractions, 
in Table || we collect our predictions for the lon- 
gitudinal polarization fractions ji, for the modes 
with two J = 1 mesons in the final state. There 
is agreement with experiment for B s — > J/ilxfri 
the only mode for which data on f L are available. Table 5: Longitudinal polarization frac- 

This is at odds with the case of a few suppressed tion h ( x10 ^ for B " dcca y s to two J = 1 

mesons 

B decays to two light vector mesons, in which the 

experimental datum is not reproduced assuming factorization. Actually, it should be no- 
ticed that naive and generalized factorization provide the same result in the case of the 
polarization fractions since they differ only for the value of which cancels in the ratio 
defining a polarization fraction. In order to modify the prediction for //_, one should either 
consider approaches in which the three polarization fractions (the longitudinal and the two 
transverse ones) are weighted by different Wilson coefficients or invoke again other mech- 
anisms such as rescattering. The first case is realized in QCD factorization and in pQCD. 
As for rescattering, it has been proposed as a solution to the puzzle of the polarization 
fractions in the case of B decays to two light vector mesons, when the considered process 
is suppressed as in the case of the penguin induced mode B <pK* [[33|] . 



Channel 


Theory 


Experiment 




51.3 ±5.8 


54.1 ± 1.7 




41.0 ±3.7 




Xd4> 


43.9 ± 4.4 





2 The quark content of fo is not completely known. Under the qq assignment, this meson might be a 
mixture of the isosinglet nn and ss (n = u, d) components. The mixing angle can be fixed using experimental 
information on, for instance, the decays J/ip —¥ 4> fo and J/ip — > ^ fo- B(J/ip — > <f>fo) = (3.2 ± 0.9) x 10 -4 , 
B(J/i/> w/ ) = (1.4 ± 0.5) x 10" 4 , which might signal a nonstrange component of fo and the consequent 
reduction of B(B S —¥ J/ipfo) by about 30% 
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3. New physics in nonleptonic B s decays: general analysis 



As mentioned in the Introduction, hints of deviations from SM predictions have recently 
been found in B s phenomenology, hence it is worth considering the effects of new physics 
in the B s sector, which may show up in mixing and/or in decay amplitudes. 

New physics in B s — B s mixing can modify the mixing phase ft s . We refer to this phase 
as to fit , which contains SM as well as NP contributions: /0f = /5f M + (3^ p . This effect 
is the same for all decay modes, and simply shifts the value of (3 S . On the other hand, 
NP in the decay amplitudes can affect various channels in different ways, even for modes 
induced by the same quark transition, as we specify in the following. 

Let us discuss the possibility that experimental results for nonleptonic B s decays de- 
viate from the predictions given in the previous Section. Such predictions rely on SU (3) 
flavour symmetry and on experimental data on corresponding Bj decays, in which no NP 
effects have been detected at the present level of accuracy. Deviations in B s decay rates 
with respect to the predictions could be due to a violation of SU(3)f symmetry, which 
is generally expected at a few percent level. However, there is the more exciting possibil- 
ity of deviations due to NP effects with small contributions in Bj oscillations and decays 
and detectable contributions in B s , an eventuality which is interesting to consider for the 
modes studied in this paper. Such modes receive contribution both from tree level and 
loop diagrams, so that one would expect NP to affect them negligibly. However, there are 
scenarios in which the contribution of new particles in loop diagrams can be competitive 
with the SM tree level diagrams. This is the case, for example, of super symmetric scenarios 
in which one loop gluino exchanges for b — > s transition could give a sizeable contribution 
to the b — > ccs induced modes. This would affect the branching ratios of such modes, 
and the CP asymmetries, since new phases could arise through the soft supersymmetry 
breaking terms. To avoid constraints on such phases from existing limits on dipole electric 
moments, one should consider flavour dependent phases. 

Here we do not focus on a specific NP model, rather we parameterize the effects of 
new physics in a general way, i.e. in terms of an amplitude, a weak and a strong phase, 
and discuss how these quantities can be constrained by experimental data on the modes 
considered above. 

In a customary notation, Af is the amplitude for B s — >■ / decay to a generic final state 
/ (CP eigenstate, common to B s and B s ) which, in our case, is of the kind M cS L 3 . The 
corresponding B s decay amplitude is denoted as Af. Being interested in CP asymmetries, 
we introduce the quantity 



in terms of which one can write the mixing induced CP asymmetry Sf and the direct CP 
asymmetry Cf. 




(3.1) 




23(A/) 
1 + |A/| : 



2 ' 




(3.2) 



3 In the case of two vectors in the final state / denotes one of the final state components being CP 
eigenstate. 
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Assuming that there is a single dominant NP amplitude (or that all NP amplitudes have 
the same weak and strong phases relative to the SM), we write: 



A f = \A s f M \ + \Af p \ e idNP e i5NP 

= \A s f M \(l + Re WNP e i5NP ) , (3.3) 

\A» P \ 

where R = is the ratio of the modulus of the NP amplitude and that of the SM one, 

while 0np(5np) is the strong (weak) NP phase with respect to the SM part. Our working 
hypothesis is that no NP affects B decays: actually, the new weak phase 5np would be the 
same in B and B s decays, depending only on the underlying quark transition, however R 
and Onp depend on the matrix elements of the operators between initial and final states, 
which can be different. 

Using the definition in eq. ( pD| ) and considering that: Af = \A s f M \{l+Re i0NP e" i5jvp ), 
we get for the CP-averaged branching fraction 

jgexp = jgSM ^ + 2R cos (^ p ) cos (S Np ) + ij2j ^ (3_ 4 ) 

where from now on we shall omit the label / in the quantities in eq. ([O]) though they 
are channel dependent (except for 5np which is the same for all the modes induced by the 
same underlying quark transition). From eq. (|3.1|), we obtain: 

1 + Re l ( NP+ s NP) y ' 



so that 



sin(2/3| // ) + 2i?cos NP sin(2/3 s e// + 5 NP ) + R 2 sin(2/3 s e// + 28 NP ) 



Sf Vf 1 + 2_Rcos Opfp cos Sjsfp + R? ' (3 ' 6) 

_ 2R sin 6 NP sin 5 NP 
f ~ 1 + 2Rcos6 NP cos 5 NP + R 2 ' ( ' ' 

r\f being the CP eigenvalue of the final state /. In absence of NP we recover the SM results 
Xf = e~ 2l/3s , Sf = — rjf sin(2/3 s ), Cf = 0. These results would also hold if NP contributes 
only to B s — B s mixing with f} a — > (3% . 



The three equations (3.4), (3.6), ( |3.7| ) allow to determine the NP parameters R, O^p, 
8np, once experimental data on B(B S — > /), Sf and Cf are available for a given final state 
/. Assuming R 1, we obtain: 

/ n,\ 

(3.8) 
(3.9) 

■» U, , M- , ■ (3 - 10) 

2 COS(PArp) COS(OArp) 




- 11 - 



where £ and Sf parametrize deviations from the SM: 



Qexp 
fiSM 



(3.11) 



-r]fSf - sin(2/3, 



cos (2/3, 



(3.12) 



In fig. HI we plot the direct CP asymmetry C / versus the mixing- induced CP asymmetry 
for several values of the strong phase 9np, using the relation ([3.8]). Once for a given channel 
/ there will be data available for (Sf, Cf) one could find a range for 6]\rp. The two panels 
in fig. H are obtained assuming different values for the B s — B s mixing phase, which in left 
panel is fixed to 2fi e J^ = 0.77 ± 0.37 rad (one of the values obtained by HFAG averaging 
the experimental results provided by the Tevatron Collaborations CDF and DO) while in 
the right panel is fixed to the SM value /3 S = 0.017 rad (since no errors are attached to the 
SM value, the various regions shrink to lines). In these figures only positive values of O^p 
have been considered, since Cf(—6) = Cf ( — \- 6 J = —Cf(8). 



C f o.o 




1.0 



0.5 



Cf o.o 



-0.5 



-1.0 







Y 


#np=- 3 


#NP- 




5 






O 




#NP=- 6 



-1.0 -0.5 0.0 0.5 1.0 

Sf 



Figure 1: Direct CP asymmetry Cf versus the mixing-induced CP asymmetry Sf for several values 
of the strong phase 9np- 



Eq. ( |3.9D allows to determine the weak phase 5np using the measured B(B S — > f) and 
Sf. To appreciate how this could be obtained, we consider the final state J/ipV^ with rjf = 1, 
and plot in Fig 5np versus Sj/^ v for S corresponding to the values B exp = 3.32 ± 1.02 
(see Table |) and B SM = 4.25 ± 0.28 (the average of the (CDSS) and (BZ) predictions in 
Table |). 

Last, we consider eq. ( 3,1C| ). In this case, all the three observables B(B S — > /), Sf and 
Cf are required to constrain R. 

As a final remark, once experimental information is available about the three observ- 
ables B(B S — > /), Sf and Cf in at least two decay modes , we could constrain also j3 s , 
since only the NP parameters R and 6np are channel-specific, while 5m p is the same for all 
the modes induced by the weak transition b — > ccs. Hence, measuring in two channels the 
six observables (B(B S — > f\), Sf x , Cf t ) and (B(B S — > $2), Sf 2 , Cf 2 ) it would be possible to 
determine Ri, 0/vpi, R2, &np,2, $np and j3 s . 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



Figure 2: Weak phase 5np versus Sj/^ in the case Onp = and for £ obtained from the values 
of B exp (left) and B SM (right) in Table 3. 



4. Conclusions 

Recent results in the B s sector require efforts to identify the most promising ways to 
unveal new physics. We have considered decay channels induced by the b — > ccs transition, 
using generalized factorization together with SU (3) f symmetry to predict their branching 
fractions in the Standard Model. Modes with a charmonium state plus rj, rj' , /o(980) 
are interesting, since they are CP eigenstates and do not require angular analyses. In 
particular, the case of /o is particularly suitable in view of its easier reconstruction in the 
7r + 7r~ mode. 

If NP affects the B s sector, it can either contribute to the AB = 2 induced mixing 
amplitude, in a channel independent way, modifying the value of the mixing phase with 
respect to its SM value, or modify the decay amplitudes, in a way that can vary from one 
channel to the other. In this case, the mixing induced CP asymmetry Sf would no more be 
equal to — rjf sin2/3 s and the direct CP asymmetry Cf would differ from zero, conditions 
which instead hold in the SM or in the case that NP affects only the oscillation process. 
Both cases can also happen. We have considered a general parametrization of new physics 
in which the second condition holds, fixing f3 s either to the average of CDF and DO results 
or to its SM value. We have shown that, with data on the branching ratio and the two 
CP asymmetries Sj and Cf for a given final state /, it will be possible to constrain NP 
parameters. 
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